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Abstract 

Using computer simulations, we studied the diffusion and structural relaxation in equilibrium 
smectic liquid crystal bulk phases of parallel hard spherocylinders. These systems exhibit a non- 
Gaussian layer-to-layer diffusion due to the presence of periodic barriers and transient cages, and 
show remarkable similarities with the behavior of out- of- equilibrium supercooled liquids. We detect 
a very slow interdayer relaxation dynamics over the whole density range of the stable smectic phase 
which spans a time interval of four time decades. The intrinsic nature of the layered structure 
yields a hopping- type diffusion which becomes more heterogeneous for higher packing fractions. In 
contrast, the in- layer dynamics is typical of a dense fluid with a relatively fast decay. Our results 
on the dynamic behavior agree well with that observed in systems of freely rotating hard rods, but 
differ quantitavely, as the height of the periodic barriers reduces to zero at the nematic-smectic 
transition for aligned rods, while it remains finite for freely rotating rods. 
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I. I. INTRODUCTION 



Non-Gaussian diffusion and dynamical heterogeneities have been shown to slow down or 
completely arrest the structural relaxation of systems, such as supercooled liquids [U El E] 
and gels jH [5]. The heterogeneous dynamics of these systems, in which individual particles 
are trapped in transient cages by neighboring particles, explains the non-exponential relax- 
ation and non-Gaussian diffusive behavior. Two scenarios are usually proposed to explain 
the non-exponential relaxation behavior: a heterogeneous scenario in which the particles 
relax exponentially at different relaxation rates, and a homogeneous scenario where the 
particles relax non-exponent ially at nearly identical rates [BJ. Recently, the analysis of dy- 
namical heterogeneities has been extended to other complex systems, such as granular media 
[7] and liquid crystals (LCs) in confined nanopores [HI IS]- Interestingly, the dynamics of LCs 
has been shown to share many features with supercooled liquids, especially in the so-called 
LC isotropic phase JTOj [11] and in the smectic LC phase [12]. The LC isotropic phase is 
a macroscopically homogeneous liquid phase, which exhibits nematic ordered domains near 
the isotropic-nematic (IN) transition. The size of these nematic domains, which are driven 
by a precursor of the nematic phase, increases upon approaching the IN-transition [13J. 
The caging behavior, i.e., the temporary localization of individual particles, as observed in 
supercooled liquids, is caused by the transient structural inhomogeneities of the nematic 
domains in the isotropic LC phase. Hence, the heterogeneous dynamics of such an isotropic 
phase with pre-nematic-order resembles that of supercooled liquids [ID] . The non-Gaussian 
behavior becomes even stronger by confining the LC in a nanoporous material [HI E] . The 
reason is that confinement modifies the dynamics of the LC particles, which yields a non- 
uniform relaxation depending on e.g., the distance to the pore surface and the pore size, 
which is indeed observed experimentally by dielectric spectroscopy [H] and quasi-elastic 
neutron scattering [HI US] • 

For the bulk smectic LCs, recent developments of experimental techniques (e.g. NMR 
coupled to strong magnetic field gradients [16], or fluorescent labeling of rods p2]), allowed 
for direct observations of non-Gaussian quasi-quantized layer-to-layer diffusion. In the light 
of these advances, Bier et al. proposed a dynamic density functional approach to study 
the self-diffusion in colloidal dispersions of infinitely elongated particles |18j . In particular, 
they investigated the effect of the local fluid structure, which temporarily cages individual 
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particles and competes with the one-dimensional "permanent" barriers due to the smectic 
layered structure. This theoretical work shows qualitative agreement with recent experi- 
ments on the self-diffusion of filamentous bacteriophage fd viruses through smectic layers 
[T7] . In both studies, the self-part of the van Hove correlation function showed clear evidence 
of an inter-layer diffusion (or permeation) occurring by discontinuous jumps of nearly one 
rod length. Simulations on freely rotating hard rods confirmed the important role of tempo- 
rary cages and permanent barriers on the non-Gaussian permeation through smectic layers, 
and revealed new insights on the relaxation behavior and cooperative motion of stringlike 
clusters [12] . 

In this work, we investigate the diffusion in bulk smectic liquid crystals of perfectly 
aligned hard spherocylinders. The phase diagram of parallel hard spherocylinders exhibit 
nematic-to-smectic (N-Sm) and smectic-to-crystal (Sm-K) phase transitions in a broad range 
of length-to-diameter ratios of the rods [HI [20]. Using computer simulations, we are able to 
study the long-time relaxation decay which approximately spans up to four time decades in 
the whole density range from the N-Sm to the Sm-K transition. We investigate the substan- 
tial differences between the fluid-like in-layer and the hopping-type inter-layer dynamics, 
which are caused by the temporary cages and permanent barriers. In addition, we study the 
heterogeneous character of the layer-to-layer dynamics, which is not exclusively a feature of 
confined LCs, and we determine the crossover between the cage regime, with the particles 
rattling around their center of mass, and the long-time diffusive regime as a function of 
density. 

II. II. MODEL AND SIMULATIONS 

Our system is composed of N = 2100 parallel hard spherocylinders with aspect ratio 
L* = L/D = 5, where D is the diameter of the hemispherical caps joined together by 
a cylindrical part of length L. Hence, the overall length of the rods is L + D. The phase 
diagram of this system displays stable nematic, smectic, and crystal phases [20]. For L* = 5, 
the smectic phase melts into a nematic phase for P* < 2.1, and crystallizes for P* > 6.3, 
where P* = Pvo/ksT is the reduced pressure with k% Boltzmann's constant and v = 
tt(D 3 /6+LD 2 /4) the molecular volume. We studied the dynamics of this system at pressures 
P* = 2.0, 2.5, 3.0, 4.0, and 5.0 corresponding to packing fractions rj = Nvq/V = 0.394, 
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0.441, 0.475, 0.525, and 0.563, respectively. We performed Monte Carlo (MC) simulation in 
a rectangular box of volume V with 7 smectic layers of aligned rods and we employ periodic 
boundary conditions. 

We first performed equilibration runs using MC simulations in the isobaric-isothermal 
(NPT) ensemble. The initial configurations for the smectic phase were obtained by expand- 
ing a solid phase at pressures P* = 2.0, 2.5, 3.0, and 4.0. The smectic phase at P* = 5.0 
was obtained by a compression run from P* = 4.0. We only performed translational moves, 
which were accepted if no overlap was detected, whereas rotational moves were not allowed. 
Volume changes have been attempted every N MC cycles by randomly changing the three 
box lengths independently. The systems were considered to be equilibrated when the packing 
fraction had reached a constant value. In the production runs, we carried out MC simula- 
tions in the canonical (NVT) ensemble, i.e., we kept the volume constant, as the collective 
moves associated with the volume changes would not correspond with Brownian dynamics. 
The maximum displacement of the MC moves was chosen in such a way that we achieve 
(i) a reasonable time of simulation, (ii) a satisfactory acceptance rate, and (ra) a suitable 
description of the Brownian motion of the particles in a colloidal suspension. To this end, 
we monitored the mean square displacement in the z- and redirections for several values of 
the maximum step size 5 max , with S maXjZ = 25 max>xy due to the anisotropy of the short-time 
self-diffusion coefficient of the rods [21]. 5 maX;Xy — D/10 and S maXjZ = D/5 were found to 
be the optimal values satisfying the above requirements. We have neglected the hydrody- 
namic effects, as it was shown recently by computer simulations of highly concentrated rod 
suspensions, that the dynamics is dominated primarily by the excluded volume and steric 
effects [22]. 

As unit of time, we have chosen r = D 2 /D tr , where D tr is the translational short-time 
diffusion coefficient, which is the isotropic average of the diffusion coefficients in the three 
space dimensions: D tr = (Dy + 2D±)/3. 

III. III. RESULTS 

The intrinsic nature of the smectic phases can be appreciated by computing the in-layer 
9±{x,y), and inter-layer g\\{z) pair correlation functions as a function of r± = \J x 2 + y 2 and 
z, respectively. In Fig. 1, we display g±(x,y) and g\\(z) for varying pressures. We find fea- 
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tures of fluid-like behavior for g±(x, y): the first peak is located at a distance approximately 
equal to one diameter length, and an exponential decay of the oscillations to one at long dis- 
tances. By contrast, g\\{z) shows pronounced periodic correlations which reveal the layered 
structure along the nematic director n. The location of the peaks corresponds to the center 
of the smectic layers, where the particle density is maximal. In Fig. 2, we show the top and 
side view of two typical configurations of the smectic LC phase at P* = 2.5 and 5.0, where 
one can clearly see the periodic structure of the smectic layers and the two-dimensional 
fluid-like structure within each layer. 




FIG. 1: (a) Inter-layer g\\(z), and (b) in-layer g±(x,y) pair correlation functions as a function of z 
and rj_ = \J x 2 + y 2 , respectively, for varying pressures as labeled. 

Upon increasing the pressure, the peaks become more pronounced indicating that the 
smectic layers are more defined and the layer spacing reduces significantly from h = 7.03/} 
at P* = 2.0 to h = QA8D at P* = 5.0. At P* = 2.0, very close to the N-Sm transition, it is 
difficult to distinguish the individual layers as the amplitude of the periodic structure reduces 
continuously, and g\\ (z) approaches a nearly flat profile in the proximity of the continuous 
N-Sm transition. At pressures P* > 4.0, the peaks become more pronounced, while the 
minima, corresponding to the interlayer spacings tends to zero. Hence, it becomes more 
difficult for the particles to diffuse from layer to layer at these high pressures. 
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FIG. 2: Top and side views of two bulk smectic liquid crystal phases observed at (a) P* = 2.5 and 
(b) P* = 5.0. 

In addition, we measure the (relative) probability tt(z) of finding a particle at position 
z with the z-axis chosen parallel to the nematic director n. We estimated the potential 
energy barrier from the Boltzmann factor tt(z) oc exp(— U{z) /ksT) as in Ref. [T7]. U(z) 
denotes the effective potential and quantifies the potential energy barrier for the layer-to- 
layer diffusion. In Fig. 3, we give U(z) at different pressures with the fitting function 
U{z) = Yli=i Ui[sm(iiz/h)} 21 , where Uo = J2i=i Ui is the potential barrier height and h the 
inter-layer spacing. We find U = 0.2k B T at P* = 2.0 and U Q = 8.3k B T at P* = 5.0. The 
potential barrier for a rod to diffuse from layer to layer increases with increasing packing 
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fraction, and becomes less steep when the N-Sm phase transition is approached at lower 
densities. 




0.5 1 1.5 2 2.5 

z/h 

FIG. 3: Effective potential U{z) in the bulk smectic phase at P* = 2.0 (O), 2.5 (■), 3.0 (O), 4.0 
(A), and 5.0 (V). The solid lines are fits. 

These results should be compared with those of Ref. [23] and the barrier heights estimated 
recently for smectic phases of freely- rotating hard rods with L* = 5 [12] . Interestingly, we 
find that the barrier height in smectic LC phases is higher for aligned hard rods than for 
freely rotating rods at the same packing fraction, i.e., we find Uq = 5.1ksT at rj = 0.508 
(P* = 3.7) for aligned rods, which should be compared with Uq = 3.5/csT for freely rotating 
rods, and Uq = 7.8/^T at r\ — 0.557 (P* = 4.8) for parallel rods to be compared with 
Uq = 7.5/c#T. The difference in barrier height decreases with increasing pressure, and we 
expect that it tends to zero upon approaching the Sm-K transition, where the freely rotating 
rods become more and more aligned. As the potential barriers are significantly lower for 
freely rotating rods at low packing fractions, we conclude that rotational degrees of freedom 
facilitate significantly the layer-to-layer diffusion. 

Barrier-free diffusion pathways can be observed even at high densities as a result of 
screw dislocations. Such structural defects create helical connections between neighboring 
smectic layers where the rods diffuse as in the nematic phase [21] . It is not an easy task to 
model screw dislocations in computer simulations, as these are not compatible with periodic 
boundary conditions. Slip boundary conditions may overcome this complication |25j . 

To study the non-Gaussian behavior of the layer-to-layer diffusion, we computed the 
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FIG. 4: Non-Gaussian parameter a.2 z (t) and mean square displacement in the direction parallel to 
the nematic director, calculated at the pressures indicated in the top frame. 

following non-Gaussian parameter [26J: 

where Az(t) = z(to + 1) — z(to) is the displacement of the rods in the z-direction in the 
time interval starting at to and ending at to + t, and (...) denotes an ensemble average 
over all particles and initial time to- For the in-layer diffusion, a similar non-Gaussian 
parameter, a2,xy{t), can be computed. In Fig. 4, we show ot2 iZ {t) at different pressures. This 
parameter quantifies the deviation from the Gaussian behavior of the probability density 
function for single-particle diffusion. At short times, the particles are freely diffusing as the 
trapping cage formed by the surrounding neighbors, is not reached yet, and hence, «2, 2 is 
basically zero. At long times, the system enters the diffusive regime and the non-Gaussian 
parameter tends to zero. At intermediate times, the motion of each particle is hampered 
by its neighbors and becomes sub-diffusive. In this time interval, a2 i2 is non-vanishing, 
indicating the development of dynamical heterogeneities. Additionally, we find that the 
peak height of a2, z (t*) at t* increases and moves to larger values of t* upon increasing the 
packing fraction, as the time for the particles to escape out of their cage increases. We also 
plot the mean square displacement (MSD) (z 2 (t)) in Fig. 4 for the same state points. The 
MSD show a clear cage-trapping plateau, which becomes more pronounced upon increasing 
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pressure. The time t* at which ct2, z displays a maximum corresponds to the end of the 
plateau observed in the MSD as can be observed clearly in Fig. 4. 

Interestingly, non-Gaussian dynamics due to cage-escape processes have also been ob- 
served in single-particle diffusion in periodic external potentials [27], 2D liquids |28j . cluster 
crystals [29], but also in glasses [301 EIJ [321 E3]- In particular, it was shown that the time 
to escape out of a cage increases because cage rearrangement (or recaging) involves a larger 
number of particles [301 Ell Ell E3] upon approaching the glass transition. at2,xy{t) ( n °t 
shown here) does not deviate significantly from zero. At the highest pressure studied, the 
peak is lower than 0.2, confirming that the in-layer relaxation dynamics is (nearly) diffusive. 
The behavior of the non-Gaussian parameter a.2, z is consistent with the theoretical predic- 
tions in systems of infinitely long parallel hard rods [H], and with the simulation results of 
freely rotating hard rods [12], although the latter exhibit higher peaks especially close to 
the smectic-to-crystal transition. Stronger deviations have been observed in colloidal sys- 
tems with short-range attractions when approaching the gel transition [31] , or in permanent 
gels where static heterogeneities give rise to a plateau at long times (i.e., the non-Gaussian 
parameter does not vanish) [HEH]. Glass transitions, by contrast, are usually characterized 
by weaker deviations [311 ES] . 

The periodic shape of the effective potential implies a hopping-type diffusion in the direc- 
tion of n, with the rods rattling around in a given layer until they overcome the free-energy 
barrier shown in Fig. 3 and jump to a neighboring layer. We quantify this layer-to-layer dif- 
fusion by calculating the self part of the Van Hove correlation function (VHF) ]36j, defined 
as: 



where 5 is the Dirac delta function. G s (z, t) measures the probability distribution for the 
z- displacements of the rods at time t + 1, given their z-positions at t . In Fig. 5, the self 
VHF is presented as a function of z at P* = 2.0 to 5.0. At P* = 3.0 to 5.0, we observe 
the appearance of peaks at distances that correspond to the center of the smectic layers in 
the z direction. No peaks are observed at P* = 2.0, where the barrier height (0.2ksT) of 
the effective potential for the layer-to-layer diffusion is sufficiently small, that there is no 
hopping-type diffusion between neighboring layers. As a result, the discontinuous diffusion 
observed in strong smectic phases, with the particles occupying quasi-discretized positions, 
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is substituted by a quasi- continuous Gaussian diffusion for weak smectic phases. At long 
times, less and less particles are still at their original positions, and consequently, the profiles 
of the VHF become almost constant in a nematic phase or periodically peaked in a smectic 
phase. The presence of peaks at distances of neighboring layers for small to intermediate 
pressures show that a significant number of particles have displaced a long distance even at 
short times. These fast-moving particles contribute to the heterogeneous dynamics of the 
system and affect its structural relaxation. In recent simulations of smectic phases of freely 
rotating hard rods, it was shown that the fast-moving particles form stringlike clusters which 
exhibit cooperative motion [T2]. We believe that this result is still valid if the orientational 
degrees of freedom are frozen out, but further investigation is needed to address this point 
in more detail. 




z/h zfli 

FIG. 5: Self part of the Van Hove function G s (z, t) at t/r = 10 (dotted lines), t/r = 100 (dashed 
lines), and t/r = 1000 (solid lines) for pressure P* = 2.0 (a), 3.0 (b), 4.0 (c), and 5.0 (d). 

Finally, we also study the structural relaxation by calculating the self-part of the inter- 
mediate scattering function: 

F s {t) =< exp[zqAr(t)] > (3) 

at wave vectors q.D = (0,0, q z ) and (q x ,Qy, 0), with q z ~ 1 and (q% + g^) 1//2 — 6, which 
correspond to the main peaks of the static structure factor. Ar(t) is the displacement of a 
particle in the time interval t. Results at pressures P* = 2.0, 3.0, and 4.0 are shown in Fig. 
6. For the pressure range P* = 2.5 — 5.0, the inter-layer dynamics shows a significantly 
slow relaxation which is characterized by a two-step decay. In the first step, which is 
relatively fast, each rod rattles around its original position without feeling the presence 
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of the surrounding neighbors. We detect an exponential decay of F s>z (t) towards a plateau 
whose height and time extension increase with pressure, as also observed for colloidal glasses 
[3]- The plateau establishes the beginning of the cage regime, where the particles start to 
interact with their nearest neighbors which form a temporary cage for the particle. The 
following step takes place at much longer times (note the logarithmic scale of Fig. 6) and 
marks the escape from the cage regime. The second decay of F s>z is well fitted by a stretched 
exponential function of the form exp[(t/tr)P], where (3 ~ 0.6 and t r is the time at which 
the intermediate scattering function decays to e _1 . The stretched exponential form of the 
relaxation decay at long times confirms the heterogeneous nature of the inter layer dynamics. 
In Fig. 7, the relaxation time is given as a function of the packing fraction and fitted with 
a power law covering 4 time decades in a density interval spanning from the N-Sm to the 
Sm-K transition. At P* = 2.0, we still observe an initial exponential decay at short times 
and a stretched exponential decay at long times, but it is very hard to detect a plateau (if 
any) due to the weak smectic character of the bulk phase and, hence, to the weak permanent 
background barriers (see Fig. 3). Nevertheless, there might be a discontinuity between two 
separated relaxation regimes, since F StZ cannot be fitted by a single scaling law. 




t/1 

FIG. 6: Self-intermediate scattering function F s (t) for the structural relaxation in the z (top) and 
xy (bottom) direction at P* = 2.0 (O), 2.5 (■), 3.0 (O), 4.0 (A), and 5.0 (v). The solid lines are 
fits. 
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By contrast, the in-layer relaxation is very fast and occurs in a single step. As can be 
observed in Fig. 7, the relaxation time, which covers basically half a time decade, does 
not change significantly as a function of density. Interestingly, the long-time decay of F S)Xy , 
which is exponential at short times, is also fitted by a stretched exponential function with 
(3 ~ 0.7 at long times. This behavior is unusual for simple fluids, where a single exponential 
decay is observed, and can also not be associated with that of a supercooled liquid as the 
characteristic cage-trapping plateau is absent. However, the in-layer dynamics in smectic 
liquid crystals is similar to that of low-density supercooled liquids which do not display a 
plateau, but only a stretched exponential relaxation at long times [57] . 
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FIG. 7: Relaxation time t r /r as a function of packing fraction r]. The solid and empty circles refer 
to the inter-layer and in-layer relaxation, respectively. The solid lines are power-law fits. 



IV. IV. CONCLUSIONS 

In conclusion, we have studied using computer simulations the diffusion of perfectly 
aligned rodlike particles in smectic liquid crystal phases, where the inter-layer dynamics 
exhibits a non-Gaussian rattling-and-jumping-type diffusion due to the simultaneous pres- 
ence of temporary cages and permanent barriers. The caging due to the mutual trapping 
of neighboring particles, slows down the diffusion, and the corresponding relaxation time 
depends strongly on the packing of the system. The presence of periodic permanent bar- 
riers with a height that increases with increasing density, are intrinsically associated to 
the layered structure of the smectic phase and determine the hopping-type layer-to-layer 
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diffusion. As detected in out-of-equilibrium colloidal systems, such as supercooled liquids, 
we found that in a given time interval some particles displace longer distances than the 
average, giving rise to a remarkable heterogeneous dynamics which results in significant 
deviations from Gaussian behavior. We quantified these heterogeneities by computing (1) 
the non-Gaussian parameter ct2,z(t), which significantly deviates from zero at high volume 
fractions; (2) the self-part of the van Hove functions, whose long tails give clear evidence 
of the presence of fast particles even at short times; (3) the mean square displacement, 
which exhibits a plateau quantifying the average life-time of the transient cages; and (4) the 
self-part of the intermediate scattering function, which characterizes the relaxation decay 
of the system. In the complete range of stability of the smectic phase, we observed a very 
slow inter-layer structural relaxation which spans over four time decades from the N-Sm to 
the Sm-K transition. We note that caging as well as the permanent barrier both contribute 
to the non-Gaussian dynamics [27]. The in-layer relaxation dynamics is very fast, but does 
not show an exponential decay at long-times, as it would be expected for simple liquids. 
The observed stretched-exponential decay corresponds to that of a low-density supercooled 
liquid where the cage effect is not sufficiently strong to yield a cage-trapping plateau. 

We observed qualitative agreement with the dynamics of freely- rotating hard rods [T2] . 
but we do find some quantitative deviations in the non-Gaussian behavior, mainly caused 
by the huge differences in the height of the potential energy barriers. To be specific, the 
height of the barrier for aligned rods tends to zero at the continuous N-Sm transition, while 
it remains finite in the case of the first order N-Sm transition of freely rotating rods, which 
changes the dynamics significantly. 
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